Efficacy of pembrolizumab and biomarker analysis in patients with WGS-based intermediate to high tumor mutational load: results from the Drug Rediscovery Protocol.

PURPOSE: To evaluate efficacy of pembrolizumab across multiple cancer types harboring different levels of Whole-Genome Sequencing (WGS)-based tumor mutational load (TML; total of non-synonymous mutations across the genome) in patients included in the Drug Rediscovery Protocol (NCT02925234). PATIENTS AND METHODS: Patients with solid, treatment-refractory, microsatellite-stable tumors were enrolled in cohort A: breast cancer TML 140-290, cohort B: tumor-agnostic cohort TML 140-290, and cohort C: tumor-agnostic cohort TML >290. Patients received pembrolizumab 200 mg every three weeks. Primary endpoint was clinical benefit (CB: objective response or stable disease (SD) ≥16 weeks). Pre-treatment tumor biopsies were obtained for WGS and RNA-sequencing. RESULTS: Seventy-two evaluable patients with 26 different histotypes were enrolled. CB rate was 13% in cohort A (3/24 with partial response (PR)), 21% in cohort B (3/24 with SD, 2/24 with PR), and 42% in cohort C (4/24 with SD, 6/24 with PR). In cohort C, neoantigen burden estimates and expression of inflammation and innate immune biomarkers were significantly associated with CB. Similar associations were not identified in cohort A and B. In cohort A, CB was significantly associated with mutations in the chromatin remodeling gene PBRM1, while in cohort B, CB was significantly associated with expression of MICA/MICB and butyrophilins. CB and clonal TML were not significantly associated. CONCLUSION: While in cohort A pembrolizumab lacked activity, cohort B and cohort C met the study's primary endpoint. Further research is warranted to refine selection of patients with tumors harboring lower TMLs and may benefit from a focus on innate immunity.

Spatial relationships in the urothelial and head and neck tumor microenvironment predict response to combination immune checkpoint inhibitors.

Immune checkpoint inhibitors (ICI) can achieve remarkable responses in urothelial cancer (UC), which may depend on tumor microenvironment (TME) characteristics. However, the relationship between the TME, usually characterized by immune cell density, and response to ICI is unclear. Here, we quantify the TME immune cell densities and spatial relationships (SRs) of 24 baseline UC samples, obtained before pre-operative combination ICI treatment, using multiplex immunofluorescence. We describe SRs by approximating the first nearest-neighbor distance distribution with a Weibull distribution and evaluate the association between TME metrics and ipilimumab+nivolumab response. Immune cell density does not discriminate between response groups. However, the Weibull SR metrics of CD8+ T cells or macrophages to their closest cancer cell positively associate with response. CD8+ T cells close to B cells are characteristic of non-response. We validate our SR response associations in a combination ICI cohort of head and neck tumors. Our data confirm that SRs, in contrast to density metrics, are strong biomarkers of response to pre-operative combination ICIs.

Evaluating the effectiveness of pre-operative diagnosis of ovarian cancer using minimally invasive liquid biopsies by combining serum human epididymis protein 4 and cell-free DNA in patients with an ovarian mass.

OBJECTIVE: To assess the feasibility of scalable, objective, and minimally invasive liquid biopsy-derived biomarkers such as cell-free DNA copy number profiles, human epididymis protein 4 (HE4), and cancer antigen 125 (CA125) for pre-operative risk assessment of early-stage ovarian cancer in a clinically representative and diagnostically challenging population and to compare the performance of these biomarkers with the Risk of Malignancy Index (RMI). METHODS: In this case-control study, we included 100 patients with an ovarian mass clinically suspected to be early-stage ovarian cancer. Of these 100 patients, 50 were confirmed to have a malignant mass (cases) and 50 had a benign mass (controls). Using WisecondorX, an algorithm used extensively in non-invasive prenatal testing, we calculated the benign-calibrated copy number profile abnormality score. This score represents how different a sample is from benign controls based on copy number profiles. We combined this score with HE4 serum concentration to separate cases and controls. RESULTS: Combining the benign-calibrated copy number profile abnormality score with HE4, we obtained a model with a significantly higher sensitivity (42% vs 0%; p<0.002) at 99% specificity as compared with the RMI that is currently employed in clinical practice. Investigating performance in subgroups, we observed especially large differences in the advanced stage and non-high-grade serous ovarian cancer groups. CONCLUSION: This study demonstrates that cell-free DNA can be successfully employed to perform pre-operative risk of malignancy assessment for ovarian masses; however, results warrant validation in a more extensive clinical study.

Combining Genomic Biomarkers to Guide Immunotherapy in Non-Small Cell Lung Cancer.

PURPOSE: The clinical value of STK11, KEAP1, and EGFR alterations for guiding immune checkpoint blockade (ICB) therapy in non-small cell lung cancer (NSCLC) remains controversial, as some patients with these proposed resistance biomarkers show durable ICB responses. More specific combinatorial biomarker approaches are urgently needed for this disease. EXPERIMENTAL DESIGN: To develop a combinatorial biomarker strategy with increased specificity for ICB unresponsiveness in NSCLC, we performed a comprehensive analysis of 254 patients with NSCLC treated with ligand programmed death-ligand 1 (PD-L1) blockade monotherapy, including a discovery cohort of 75 patients subjected to whole-genome sequencing (WGS), and an independent validation cohort of 169 patients subjected to tumor-normal large panel sequencing. The specificity of STK11/KEAP1/EGFR alterations for ICB unresponsiveness was assessed in the contexts of a low (<10 muts/Mb) or high (≥10 muts/Mb) tumor mutational burden (TMB). RESULTS: In low TMB cases, STK11/KEAP1/EGFR alterations were highly specific biomarkers for ICB resistance, with 0/15 (0.0%) and 1/34 (2.9%) biomarker-positive patients showing treatment benefit in the discovery and validation cohorts, respectively. This contrasted with high TMB cases, where 11/13 (85%) and 15/34 (44%) patients with at least one STK11/KEAP1/EGFR alteration showed durable treatment benefit in the discovery and validation cohorts, respectively. These findings were supported by analyses of progression-free survival and overall survival. CONCLUSIONS: The unexpected ICB responses in patients carrying resistance biomarkers in STK11, KEAP1, and EGFR were almost exclusively observed in patients with a high TMB. Considering these alterations in context, the TMB offered a highly specific combinatorial biomarker strategy for limiting overtreatment in NSCLC.

Distinct spatiotemporal dynamics of CD8+ T cell-derived cytokines in the tumor microenvironment.

Cells in the tumor microenvironment (TME) influence each other through secretion and sensing of soluble mediators, such as cytokines and chemokines. While signaling of interferon γ (IFNγ) and tumor necrosis factor α (TNFα) is integral to anti-tumor immune responses, our understanding of the spatiotemporal behavior of these cytokines is limited. Here, we describe a single cell transcriptome-based approach to infer which signal(s) an individual cell has received. We demonstrate that, contrary to expectations, CD8+ T cell-derived IFNγ is the dominant modifier of the TME relative to TNFα. Furthermore, we demonstrate that cell pools that show abundant IFNγ sensing are characterized by decreased expression of transforming growth factor ß (TGFß)-induced genes, consistent with IFNγ-mediated TME remodeling. Collectively, these data provide evidence that CD8+ T cell-secreted cytokines should be categorized into local and global tissue modifiers, and describe a broadly applicable approach to dissect cytokine and chemokine modulation of the TME.

MYC is a clinically significant driver of mTOR inhibitor resistance in breast cancer.

Targeting the PI3K-AKT-mTOR pathway is a promising therapeutic strategy for breast cancer treatment. However, low response rates and development of resistance to PI3K-AKT-mTOR inhibitors remain major clinical challenges. Here, we show that MYC activation drives resistance to mTOR inhibitors (mTORi) in breast cancer. Multiomic profiling of mouse invasive lobular carcinoma (ILC) tumors revealed recurrent Myc amplifications in tumors that acquired resistance to the mTORi AZD8055. MYC activation was associated with biological processes linked to mTORi response and counteracted mTORi-induced translation inhibition by promoting translation of ribosomal proteins. In vitro and in vivo induction of MYC conferred mTORi resistance in mouse and human breast cancer models. Conversely, AZD8055-resistant ILC cells depended on MYC, as demonstrated by the synergistic effects of mTORi and MYCi combination treatment. Notably, MYC status was significantly associated with poor response to everolimus therapy in metastatic breast cancer patients. Thus, MYC is a clinically relevant driver of mTORi resistance that may stratify breast cancer patients for mTOR-targeted therapies.

Multi-omics analysis reveals distinct non-reversion mechanisms of PARPi resistance in BRCA1- versus BRCA2-deficient mammary tumors.

BRCA1 and BRCA2 both function in DNA double-strand break repair by homologous recombination (HR). Due to their HR defect, BRCA1/2-deficient cancers are sensitive to poly(ADP-ribose) polymerase inhibitors (PARPis), but they eventually acquire resistance. Preclinical studies yielded several PARPi resistance mechanisms that do not involve BRCA1/2 reactivation, but their relevance in the clinic remains elusive. To investigate which BRCA1/2-independent mechanisms drive spontaneous resistance in vivo, we combine molecular profiling with functional analysis of HR of matched PARPi-naive and PARPi-resistant mouse mammary tumors harboring large intragenic deletions that prevent reactivation of BRCA1/2. We observe restoration of HR in 62% of PARPi-resistant BRCA1-deficient tumors but none in the PARPi-resistant BRCA2-deficient tumors. Moreover, we find that 53BP1 loss is the prevalent resistance mechanism in HR-proficient BRCA1-deficient tumors, whereas resistance in BRCA2-deficient tumors is mainly induced by PARG loss. Furthermore, combined multi-omics analysis identifies additional genes and pathways potentially involved in modulating PARPi response.

ecpc: an R-package for generic co-data models for high-dimensional prediction.

BACKGROUND: High-dimensional prediction considers data with more variables than samples. Generic research goals are to find the best predictor or to select variables. Results may be improved by exploiting prior information in the form of co-data, providing complementary data not on the samples, but on the variables. We consider adaptive ridge penalised generalised linear and Cox models, in which the variable-specific ridge penalties are adapted to the co-data to give a priori more weight to more important variables. The R-package ecpc originally accommodated various and possibly multiple co-data sources, including categorical co-data, i.e. groups of variables, and continuous co-data. Continuous co-data, however, were handled by adaptive discretisation, potentially inefficiently modelling and losing information. As continuous co-data such as external p values or correlations often arise in practice, more generic co-data models are needed. RESULTS: Here, we present an extension to the method and software for generic co-data models, particularly for continuous co-data. At the basis lies a classical linear regression model, regressing prior variance weights on the co-data. Co-data variables are then estimated with empirical Bayes moment estimation. After placing the estimation procedure in the classical regression framework, extension to generalised additive and shape constrained co-data models is straightforward. Besides, we show how ridge penalties may be transformed to elastic net penalties. In simulation studies we first compare various co-data models for continuous co-data from the extension to the original method. Secondly, we compare variable selection performance to other variable selection methods. The extension is faster than the original method and shows improved prediction and variable selection performance for non-linear co-data relations. Moreover, we demonstrate use of the package in several genomics examples throughout the paper. CONCLUSIONS: The R-package ecpc accommodates linear, generalised additive and shape constrained additive co-data models for the purpose of improved high-dimensional prediction and variable selection. The extended version of the package as presented here (version number 3.1.1 and higher) is available on ( https://cran.r-project.org/web/packages/ecpc/ ).

A living biobank of patient-derived ductal carcinoma in situ mouse-intraductal xenografts identifies risk factors for invasive progression.

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive breast cancer (IBC). Due to a lack of biomarkers able to distinguish high- from low-risk cases, DCIS is treated similar to early IBC even though the minority of untreated cases eventually become invasive. Here, we characterized 115 patient-derived mouse-intraductal (MIND) DCIS models reflecting the full spectrum of DCIS observed in patients. Utilizing the possibility to follow the natural progression of DCIS combined with omics and imaging data, we reveal multiple prognostic factors for high-risk DCIS including high grade, HER2 amplification, expansive 3D growth, and high burden of copy number aberrations. In addition, sequential transplantation of xenografts showed minimal phenotypic and genotypic changes over time, indicating that invasive behavior is an intrinsic phenotype of DCIS and supporting a multiclonal evolution model. Moreover, this study provides a collection of 19 distributable DCIS-MIND models spanning all molecular subtypes.

PD-L1 blockade in combination with carboplatin as immune induction in metastatic lobular breast cancer: the GELATO trial.

Invasive lobular breast cancer (ILC) is the second most common histological breast cancer subtype, but ILC-specific trials are lacking. Translational research revealed an immune-related ILC subset, and in mouse ILC models, synergy between immune checkpoint blockade and platinum was observed. In the phase II GELATO trial ( NCT03147040 ), patients with metastatic ILC were treated with weekly carboplatin (area under the curve 1.5 mg ml-1 min-1) as immune induction for 12 weeks and atezolizumab (PD-L1 blockade; triweekly) from the third week until progression. Four of 23 evaluable patients had a partial response (17%), and 2 had stable disease, resulting in a clinical benefit rate of 26%. From these six patients, four had triple-negative ILC (TN-ILC). We observed higher CD8+ T cell infiltration, immune checkpoint expression and exhausted T cells after treatment. With this GELATO trial, we show that ILC-specific clinical trials are feasible and demonstrate promising antitumor activity of atezolizumab with carboplatin, particularly for TN-ILC, and provide insights for the design of highly needed ILC-specific trials.

Codon-specific KRAS mutations predict survival benefit of trifluridine/tipiracil in metastatic colorectal cancer.

Genomics has greatly improved how patients with cancer are being treated; however, clinical-grade genomic biomarkers for chemotherapies are currently lacking. Using whole-genome analysis of 37 patients with metastatic colorectal cancer (mCRC) treated with the chemotherapy trifluridine/tipiracil (FTD/TPI), we identified KRAS codon G12 (KRASG12) mutations as a potential biomarker of resistance. Next, we collected real-world data of 960 patients with mCRC receiving FTD/TPI and validated that KRASG12 mutations were significantly associated with poor survival, also in analyses restricted to the RAS/RAF mutant subgroup. We next analyzed the data of the global, double-blind, placebo-controlled, phase 3 RECOURSE trial (n = 800 patients) and found that KRASG12 mutations (n = 279) were predictive biomarkers for reduced overall survival (OS) benefit of FTD/TPI versus placebo (unadjusted interaction P = 0.0031, adjusted interaction P = 0.015). For patients with KRASG12 mutations in the RECOURSE trial, OS was not prolonged with FTD/TPI versus placebo (n = 279; hazard ratio (HR) = 0.97; 95% confidence interval (CI) = 0.73-1.20; P = 0.85). In contrast, patients with KRASG13 mutant tumors showed significantly improved OS with FTD/TPI versus placebo (n = 60; HR = 0.29; 95% CI = 0.15-0.55; P < 0.001). In isogenic cell lines and patient-derived organoids, KRASG12 mutations were associated with increased resistance to FTD-based genotoxicity. In conclusion, these data show that KRASG12 mutations are biomarkers for reduced OS benefit of FTD/TPI treatment, with potential implications for approximately 28% of patients with mCRC under consideration for treatment with FTD/TPI. Furthermore, our data suggest that genomics-based precision medicine may be possible for a subset of chemotherapies.

Tumour-educated platelets for breast cancer detection: biological and technical insights.

BACKGROUND: Studies have shown that blood platelets contain tumour-specific mRNA profiles tumour-educated platelets (TEPs). Here, we aim to train a TEP-based breast cancer detection classifier. METHODS: Platelet mRNA was sequenced from 266 women with stage I-IV breast cancer and 212 female controls from 6 hospitals. A particle swarm optimised support vector machine (PSO-SVM) and an elastic net-based classifier (EN) were trained on 71% of the study population. Classifier performance was evaluated in the remainder (29%) of the population, followed by validation in an independent set (37 cases and 36 controls). Potential confounding was assessed in post hoc analyses. RESULTS: Both classifiers reached an area under the curve (AUC) of 0.85 upon internal validation. Reproducibility in the independent validation set was poor with an AUC of 0.55 and 0.54 for the PSO-SVM and EN classifier, respectively. Post hoc analyses indicated that 19% of the variance in gene expression was associated with hospital. Genes related to platelet activity were differentially expressed between hospitals. CONCLUSIONS: We could not validate two TEP-based breast cancer classifiers in an independent validation cohort. The TEP protocol is sensitive to within-protocol variation and revision might be necessary before TEPs can be reconsidered for breast cancer detection.

γδ T cells are effectors of immunotherapy in cancers with HLA class I defects.

DNA mismatch repair-deficient (MMR-d) cancers present an abundance of neoantigens that is thought to explain their exceptional responsiveness to immune checkpoint blockade (ICB)1,2. Here, in contrast to other cancer types3-5, we observed that 20 out of 21 (95%) MMR-d cancers with genomic inactivation of ß2-microglobulin (encoded by B2M) retained responsiveness to ICB, suggesting the involvement of immune effector cells other than CD8+ T cells in this context. We next identified a strong association between B2M inactivation and increased infiltration by γδ T cells in MMR-d cancers. These γδ T cells mainly comprised the Vδ1 and Vδ3 subsets, and expressed high levels of PD-1, other activation markers, including cytotoxic molecules, and a broad repertoire of killer-cell immunoglobulin-like receptors. In vitro, PD-1+ γδ T cells that were isolated from MMR-d colon cancers exhibited enhanced reactivity to human leukocyte antigen (HLA)-class-I-negative MMR-d colon cancer cell lines and B2M-knockout patient-derived tumour organoids compared with antigen-presentation-proficient cells. By comparing paired tumour samples from patients with MMR-d colon cancer that were obtained before and after dual PD-1 and CTLA-4 blockade, we found that immune checkpoint blockade substantially increased the frequency of γδ T cells in B2M-deficient cancers. Taken together, these data indicate that γδ T cells contribute to the response to immune checkpoint blockade in patients with HLA-class-I-negative MMR-d colon cancers, and underline the potential of γδ T cells in cancer immunotherapy.

IL-5-producing CD4+ T cells and eosinophils cooperate to enhance response to immune checkpoint blockade in breast cancer.

Immune checkpoint blockade (ICB) has heralded a new era in cancer therapy. Research into the mechanisms underlying response to ICB has predominantly focused on T cells; however, effective immune responses require tightly regulated crosstalk between innate and adaptive immune cells. Here, we combine unbiased analysis of blood and tumors from metastatic breast cancer patients treated with ICB with mechanistic studies in mouse models of breast cancer. We observe an increase in systemic and intratumoral eosinophils in patients and mice responding to ICB treatment. Mechanistically, ICB increased IL-5 production by CD4+ T cells, stimulating elevated eosinophil production from the bone marrow, leading to systemic eosinophil expansion. Additional induction of IL-33 by ICB-cisplatin combination or recombinant IL-33 promotes intratumoral eosinophil infiltration and eosinophil-dependent CD8+ T cell activation to enhance ICB response. This work demonstrates the critical role of eosinophils in ICB response and provides proof-of-principle for eosinophil engagement to enhance ICB efficacy.

Assessment of Predictive Genomic Biomarkers for Response to Cisplatin-based Neoadjuvant Chemotherapy in Bladder Cancer.

Cisplatin-based neoadjuvant chemotherapy (NAC) followed by radical cystectomy is recommended for patients with muscle-invasive bladder cancer (MIBC). It has been shown that somatic deleterious mutations in ERCC2, gain-of-function mutations in ERBB2, and alterations in ATM, RB1, and FANCC are correlated with pathological response to NAC in MIBC. The objective of this study was to validate these genomic biomarkers in pretreatment transurethral resection material from an independent retrospective cohort of 165 patients with MIBC who subsequently underwent NAC and radical surgery. Patients with ypT0/Tis/Ta/T1N0 disease after surgery were defined as responders. Somatic deleterious mutations in ERCC2 were found in nine of 68 (13%) evaluable responders and two of 95 (2%) evaluable nonresponders (p = 0.009; FDR = 0.03). No correlation was observed between response and alterations in ERBB2 or in ATM, RB1, or FANCC alone or in combination. In an exploratory analysis, no additional genomic alterations discriminated between responders and nonresponders to NAC. No further associations were identified between the aforementioned biomarkers and pathological complete response (ypT0N0) after surgery. In conclusion, we observed a positive association between deleterious mutations in ERCC2 and pathological response to NAC, but not overall survival or recurrence-free survival. Other previously reported genomic biomarkers were not validated. PATIENT SUMMARY: It is currently unknown which patients will respond to chemotherapy before definitive surgery for bladder cancer. Previous studies described several gene mutations in bladder cancer that correlated with chemotherapy response. This study confirmed that patients with bladder cancer with a mutation in the ERCC2 gene often respond to chemotherapy.

Homologous Recombination Deficiency and Cyclin E1 Amplification Are Correlated with Immune Cell Infiltration and Survival in High-Grade Serous Ovarian Cancer.

BACKGROUND: How molecular profiles are associated with tumor microenvironment (TME) in high-grade serous ovarian cancer (HGSOC) is incompletely understood. Therefore, we analyzed the TME and molecular profiles of HGSOC and assessed their associations with overall survival (OS). METHODS: Patients with advanced-stage HGSOC treated in three Dutch hospitals between 2008-2015 were included. Patient data were collected from medical records. BRCA1/2 mutation, BRCA1 promotor methylation analyses, and copy number variations were used to define molecular profiles. Immune cells were assessed with immunohistochemical staining. RESULTS: 348 patients were categorized as BRCA mutation (BRCAm) (BRCAm or promotor methylation) (30%), non-BRCA mutated HRD (19%), Cyclin E1 (CCNE1)-amplification (13%), non-BRCAmut HRD and CCNE1-amplification (double classifier) (20%), and no specific molecular profile (NSMP) (18%). BRCAm showed highest immune cell densities and CCNE1-amplification lowest. BRCAm showed the most favorable OS (52.5 months), compared to non-BRCAmut HRD (41.0 months), CCNE1-amplification (28.0 months), double classifier (27.8 months), and NSMP (35.4 months). Higher immune cell densities showed a favorable OS compared to lower, also within the profiles. CD8+, CD20+, and CD103+ cells remained associated with OS in multivariable analysis. CONCLUSIONS: Molecular profiles and TME are associated with OS. TME differs per profile, with higher immune cell densities showing a favorable OS, even within the profiles. HGSOC does not reflect one entity but comprises different entities based on molecular profiles and TME.

Extensive androgen receptor enhancer heterogeneity in primary prostate cancers underlies transcriptional diversity and metastatic potential.

Androgen receptor (AR) drives prostate cancer (PCa) development and progression. AR chromatin binding profiles are highly plastic and form recurrent programmatic changes that differentiate disease stages, subtypes and patient outcomes. While prior studies focused on concordance between patient subgroups, inter-tumor heterogeneity of AR enhancer selectivity remains unexplored. Here we report high levels of AR chromatin binding heterogeneity in human primary prostate tumors, that overlap with heterogeneity observed in healthy prostate epithelium. Such heterogeneity has functional consequences, as somatic mutations converge on commonly-shared AR sites in primary over metastatic tissues. In contrast, less-frequently shared AR sites associate strongly with AR-driven gene expression, while such heterogeneous AR enhancer usage also distinguishes patients' outcome. These findings indicate that epigenetic heterogeneity in primary disease is directly informative for risk of biochemical relapse. Cumulatively, our results illustrate a high level of AR enhancer heterogeneity in primary PCa driving differential expression and clinical impact.